The effect of a small magnetic flux on the metal-insulator transition 
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We numerically show that very small magnetic flux can significantly shift the metal-insulator 
transition point in a disordered electronic system. The shift we observe for the 3d Anderson model 
obeys a power law as predicted by Larkin and Khmel'nitskii (1981). We compute the exponent 
and find good agreement with the prediction. However, the power law holds only for much smaller 
magnetic fields than has been previously assumed, and is accompanied by a large prefactor, leading to 
a surprising strong dependence of the transition point on the applied magnetic field. Furthermore, we 
show that the critical level-spacing distribution is identical in the presence and absence of a magnetic 
field if hard-wall boundary conditions are applied. This result is surprising since both cases belong 
to different universality classes and different distributions have been reported for periodic boundary 
conditions. 

PACS numbers: 71.30.-|-h, 75.47.-m, 72.15.Rn 
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I. INTRODUCTION 

The interest in measuring small magnetic fields is 
driven by the possibility to build smaller magnetic stor- 
age devices with high capacity. Much progress in the 
understanding of the magnetic properties of condensed 
matter has been achieved in the last decades especially 
in thin layers^. Some of the results were honored by a 
Nobel price in physics in 2007 for the discovery of the 
giant magneto resistance by Fert and Griinberg. 

In this paper we want to take a different path, study- 
ing the influence of a magnetic field on a metal-insulator 
transition (Anderson transition)^ rather than on the 
magnetization. Already in the eighties, Larkin and 
Khmel'nitskii^ estimated the shift of the metal-insulator 
transition point, i.e., the critical disorder Wc for small 
magnetic fluxes as 

Wc{(l)) - Wc{0) = A(j)'^/'' with P = 1/2 (1) 

and v the scaling exponent of the localization length, 
X{W) ^ {W - Wc)-". The value z/ « 1.43 was numeri- 
cally calculated for large and random magnetic fluxes in 
3D cubic latticesi. The prediction ([T]) was verified with 
different sophisticated analytical studies^'^'^'^'^ '^^i^^ . On 
the other hand, numerical work concentrated on the ef- 
fect of large or random magnetic fluxes^. Computation 
of the shift of the mobility edge as function of the mag- 
netic field shift was performed only oncei^, with large 
error bars and for relatively large fluxes. However, since 
Eq. Il]) should be correct only for small fluxes a more 
thorough numerical analysis is needed. 

Here we show that for small magnetic fluxes the shift of 
the critical disorder Wc is a power law following Eq. ([1]) 
very closely. Deviations occur for </> > 0.03, in particu- 
lar for the (j) values previously considered^^. Moreover, 
our numerical results show a surprisingly large prefactor 
A in Eq. ([1]), which leads to a very large shift even for 
small values of (f>. Basically, half of the total shift in Wc 
takes place for fluxes which are more than hundred times 



smaller than the ones previously considered. This could 
possibly be exploited for devices detecting small mag- 
netic fields. Various experiments have demonstrated the 
possibility to reproduce properties seen in the Anderson 
model including effects of magnetic field o^'^i^'^'^^ . 

In addition to the shift of the critical disorder Wc in- 
duced by > a change in the critical level-spacing 
distribution with magnetic field was reportediS. These 
results were obtained with periodic boundary condi- 
tions. On the other hand, it is known that a change 
of boundary conditions also changes the critical level- 
spacing distributioniiii^. Here we show that for hard- 
wall boundary conditions one finds identical critical level- 
spacing distribution functions in the presence and ab- 
sence of a magnetic flux, even though the magnetic field 
changes the universality class of the system. We also ex- 
amine the second moment of the critical level spacing dis- 
tributions and confirm its dependence on the boundary 
condition^^. Furthermore, we examine its dependence on 
the magnetic fiux for periodic and hard wall boundary 
conditions and find no changes in both cases. 

The outline of the paper is as follows. In Section |TT] 
we introduce the model and also describe the finite-size 
fitting procedures we employed to determine the criti- 
cal disorders Wc- In Section Hill we discuss the main re- 
sults presenting the phase diagram and comparing with 
Eq. ^ . Section [IV] is devoted to the effect of bound- 
ary conditions on the critical level spacing distribution 
and its second moment. Section Ivl gives conclusion and 
outlook. 



II. MODEL AND FINITE-SIZE SCALING 
APPROACH 



To study metal-insulator transitions (MIT) on a 
3D simple cubic lattice we consider the tight-binding 
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Hamiltonian^, 

-ff = ^ eiojoi - ^ ij^iojoi, (2) 

where the first part represents the disordered on-site 
(node) potential (homogeneous distribution —W/2 < 
Ci < W/2) and the second part describes the transfer 
between neighboring sites («, j). The transfer probabihty 
is given by tij — exp^iipij) for neighboring lattice points 
where ipi^ is the phase accumulation of the hopping elec- 
tron. A constant magnetic field B is applied by choosing 
phases such that the flux cj) per plaquete (unit square) 
in the xy-plane is in between 8 x 10~^ and 0.25 in units 
of the flux quantum (po = hc/e. The lowest value is de- 
rived from the smallest system size L/a — 14 we consider 
here, with a the lattice constant defining the length scale. 
For lower fluxes the magnetic length Ln/a — l/(-\/27r0) 
would become smaller than L, and the electrons will not 
accumulate a 27r phase shift. The largest flux we consider 
is given by the symmetry of the problemiS. The relation 
between the flux per plaquete and a given magnetic field 
depends on the lattice constant a, (f> = eBa^/h. We 
choose hard wall boundary conditions in most calcula- 
tions to avoid discontinuities in the phase accumulation 
of the hopping electrons due to the magnetic fleld. 

To calculate the eigenvalues of 3D systems with the 
Hamiltonian ^ we use an iterative solver based on a gen- 
eral Davidson and Olsen algorithm^^'^ where the matrix 
vector multiplication is performed by Intel MKL Par- 
diso. We calculate six eigenvalues around \E\ — for 
each configuration with linear system sizes ranging from 
L/a = 14 to L/a — 40 lattice points, accumulating at 
least 2 X 10^ eigenvalues for each size. 

It is expected that a second-order phase transition 
occurs for a given critical disorder Wc- For an infi- 
nite system this transition is characterized by the di- 
vergence of the correlation length A. The divergence is 
described by A(iy) ~ {W — Wc)~'^ with a critical expo- 
nent v. We study this transition by analyzing the level- 
spacing distribution P{s) of consecutive eigenvalues Ei, 
with s = {Ei — Ei^i)/AE and the mean level spacing 
AE. 

In the limit of infinite system size the level-spacing 
distribution P{s) of a disordered system corresponds to 
a random-matrix theory result if the eigenfunctions are 
extended. For systems which obey time-reversal symme- 
try the corresponding random-matrix ensemble is the or- 
thogonal Gaussian ensembles (GOE), well approximated 
by the Wigner surmise, Pgoe(s) = (7r/2)s exp(— 7rs^/4) 
(dashed red line in Fig. [1]). A magnetic field breaks 
the time reversal symmetry, changing the ensemble to 
a unitary Gaussian ensemble (GUE), for which the level- 
spacing distribution is approximated by i-GUE(s) = 
(32/7r2)s2exp(-4sV7r) (dotted red line in Fig. [T]). In 
contrast to the metallic phase, the uncorrelated energy 
levels of localized states are characterized by a Poisson 
distribution, Pp{s) = exp(— s) (dashed dotted purple line 
in Fig.[T|), independent of the symmetry. 
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FIG. 1: (Color online) The level spacing distribution P{s) for 
a system of size L/a — 40 and magnetic flux per plaquete 
(j> = 10~^. In the metallic phase (red curves) the level spacing 
follows the GUE distribution (dotted red curve) rather than 
the GOE distribution (dashed red curve). A clear transition 
from GUE to Poisson (dash-dotted purple curve) behavior 
is observed as the disorder increases from W — 8 (continu- 
ous red curve) to = 23 (continuous purple curve). Inset: 
System size dependence of Jo = (s^)/2 for W = 10. With 
increasing system size the GUE value (dotted line) is clearly 
reached asymptotically rather than the GOE value (dashed 
line). 

For finite systems P{s) is between Pgoe (or Pgue with 
a magnetic field) and Pp(s). However, it approaches one 
of them with increasing system size, remaining system- 
size independent only at the transition point, i.e. for 
W = Wc- To determine the critical disorder we study the 
system-size (L) dependence, of the second moment of the 
level spacings, Iq = (s^)/2. From one-parameter finite- 
size scaling arguments^i^i we expect that Iq will depend 
on disorder W and lattice size L \iW ^ Wc, but become 
independent oi L aiW — Wc- The second moments of a 
Poissonian, GOE or GUE distribution can be calculated 
via s'^P{s), yielding /o,p = 1, /q.goe — 0.637 (dashed 
line in the inset of Fig.[T]) and /q.gue = 0.589 (dotted line 
in the inset of Fig. [J). 

Since we choose hard boundary conditions we need to 
take care of finite-size corrections due to irrelevant sur- 
face effects^. Therefore we fit our data to a scaling form 
including these irrelevant surface effects decaying with 
system size as a power law, Iq = F(^L^/'', SL^) with the 
critical exponent v, the relevant scaling variable ^, the 
leading irrelevant variable S and the leading irrelevant 
exponent y- For finite system sizes L no phase transition 
takes place and F is a smooth analytical function. After 
expanding F to first order one gets: 

/o(0,iy,L) = Fo(vI/Li/-) + SL^i^i(vI/Li/-). (3) 

where Fq and Fi are analytical functions and are ex- 
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panded to third order in the analysis of the numerical 
results. The relevant and irrelevant scaling variables are 
expanded in power series of the dimensionless disorder 
w = {Wc — W)/Wc- The relevant variable is expanded 
to first order, "^{w) — '^i{Wc — W)/Wc, whereas we ex- 
pand to zeroth order the irrelevant variable, S(w) — Sq. 
In total there are eleven independent fit parameters with 
Wc and 1/ being the interesting ones (see Refj2S, for the 
details of the procedure). 

Figure [T] shows the level-spacing distributions of a sys- 
tem of size L/a = 40 and magnetic flux (p = 10""^ per 
plaquete. A transition can be observed from a metal- 
lic phase, where the level-spacing distribution is close to 
GUE, for W = 8 (continuous red line), to an insulating 
phase with P{s) close to Poisson, for W = 23 (contin- 
uous purple line). When comparing the continuous red 
line with the dotted red curve, a clear similarity can be 
observed. This means that the metallic state is simi- 
lar to a random matrix GUE (Pgue(s)) and not GOE 
(Pgoe(s)). We do not observe a smooth crossover be- 
tween GOE and GUE for linear system sizes larger than 
Lh- The same asymptotic behavior can also be seen in 
the size dependence of the second moment of the level 
spacings Iq for a metallic state. For large system sizes 
and low disorder {W = 10), Iq reaches the GUE value 
asymptotically as depicted in the inset of Fig. [T] The 
critical (system-size independent) level-spacing distribu- 
tion occurs at Wc « 17 (continuous hght blue line). This 
distribution will be discussed in more detail in Section 
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Since the actual form of the critical level-spacing dis- 
tribution is not known we study Iq to extract the crit- 
ical disorder values Wc and the corresponding la.c- In 
the procedure we fit our results of Iq to Eq. ^ using 
a least-square method, i.e., a state of the art non- linear 
least-square fitting algorithm implemented in Mathemat- 
ica. To achieve better accuracy we do not fit the results 
for each flux (f> separately but simultaneously calculate 
a non-linear fit over certain ranges of fluxes. Such an 
approach is possible since Wc is the only parameter in 
Eq. ^ changing significantly with (j) for small values of 
(j). To check the stability of the fits we compare several 
fits over different ranges of fiux between = 8 x 10"'' 
and = 3 X lO'^. For large flux {(j) = 0.1 and 0.25) as 
well as for no flux (0 — 0) we perform separate fits, since 
other parameters than Wc (in particularly change for 
the GOE ensemble {4> = 0) and for large (f>. 

Figure HJa) shows results of such a non-linear fit for 
(p = 10~^ together with the simulation results for four 
system sizes {L/a = 14,20,30,40). As expected and al- 
ready indicated in the inset of Fig. [l] Iq drops with in- 
creasing system sizes reaching towards /o,gue for large 
system sizes if < Wc- For W > Wc, on the other 
hand, Iq increases with system sizes towards /o,p. At the 
critical disorder Wc, Iq is system-size independent. The 
lines are the result of the simultaneous fit. Although it 
seems that all lines cross atW^ 16.8 this is not the case 
as can be seen in the inset. This is due to the non-relevant 




FIG. 2: (Color online) The localization parameter 7o 



')/2 



versus disorder W for system sizes L/a = 14 (red circles), 
L/a — 20 (green squares), L/a = 30 (light blue diamonds), 
and L/a = 40 (blue stars) where a magnetic flux of (j) — 
10"'^ is applied. The lines correspond to fits of Eq. Q. (a) 
The symbols correspond to Jq. A transition from extended 
states for small W to localized ones for large W is seen for all 
sample sizes. Nevertheless, the critical crossing point cannot 
be observed directly due to finite size effects. Inset: the region 
around W = 16.8 zoomed in; there is no single crosspoint. (b) 
Corrected values of Iq where the influence of the irrelevant 
variables is subtracted. A clear transition can be seen at 
Wc ~ 17: all lines cross at one distinct critical disorder. 



variables EL^Fi scaling with system size L. Figure [IJb) 
shows that subtracting the non-relevant variables leads to 
a nice crossover for Wc ~ 17 indicating a system-size in- 
dependent critical value of Iqc ~ 0.83. Taking the whole 
range of fluxes it is possible to draw a phase diagram of 
fluxes and critical disorders. 
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TABLE I: Result of single fits for the interesting parameters. 









0.1 


0.25 




f6.59 


18.21 


18.30 


V 


1.54 


1.42 


1.49 


Io,c 


0.83 


0.83 


0.81 



III. PHASE DIAGRAM AND SCALING 
BEHAVIOR OF THE TRANSITION 

We have performed extensive numerical simulations for 
the entire range of magnetic fluxesi^ from (p — to 0.25. 
The data is analyzed according to the description given 
in the previous section. The results for the extreme cases, 
(j) — 0, 0.1, and 0.25, are presented in Tab.Hl For 4> = 
we retrieve the well-known GOE values for Wc and 
For = 0.1 and — 0.25 we also reproduce previous 
results^ii^. Interestingly, the critical value /q.c does not 
change upon applying the flux, which is also in agree- 
ment with previous work considering periodic boundary 
conditions^ii^iii. 

Figure [3] summarizes our results for small fluxes, 4> — 
8 X lO^'^, 10~^ 1.5 X 10-^, 2 X 10~3, 4 x 10"^ 7 x 10"^ 
10~^, and 3 x lO^'^, based on a simultaneous fit. The 
dashed line marks the flux where the magnetic length 
Lh/cl — ^2-K<p coincides with our smallest system size 
L/a = 14. The values of the most important parameters 
Wc, V and /qc are also reported in Tab. [TTl They are 
means of the results achieved from different starting val- 
ues in the fitting procedure; the errors are the absolute 
deviations from the means. The result for the critical 
exponent at small fluxes v = 1.41 is slightly lower but 
within the error bar range of a previous numerical studj>^ 
(y — 1.43 ±0.04). This result is similar to our results for 
larger fluxes (0 — 0.1 and 0.25). 

Our main goal , however, is to study the dependence 
of the critical disorder Wc on 0. Since the fitting ap- 
proach does not use as a fitting parameter, the power- 
law behavior shown in Fig. [3] confirms the prediction 
(fT|) by Larkin and Khmel'nitskii^ discussed in the intro- 
duction. Although the individual error bars for Wc((/)) 
are still significantly large, the power-law behavior for 
VFc(</') comes out naturally, even though Eq. ^ is not 
used in our scaling approach. The calculated exponent 
/3 = 0.49±0.08 is in perfect agreement with the predicted 
exponent /3 = l/2i^i^i^i^i^ i^°i^^ The error bar represents 
the deviations we get in the exponent when fitting the 
different results of Wc{<fy for the different initial fit vari- 
ables. When different regimes of are used for the si- 
multaneous fit, all tVc(0) values shift slightly in some 
direction, but /3 stays the same. Therefore /3 is more 
exact than the error bars in Fig. [3] might indicate. 

The inset of Fig. [3] shows the critical disorder Wc{(l>) 
(indicated by squares) as function of 0. The red line is the 
same power-law fit as in the main panel. Although the 




FIG. 3: (Color online) The shift of the critical disorder Wc 
as function of the applied magnetic field. The symbols corre- 
spond to the means obtained in a simultaneous non-linear fit 
of Eq. (O with fluxes between 8 x 10"" and 3 x 10"^. The 
error bars are the absolute deviation from the mean for dif- 
ferent initial fitting parameters. The red line is a linear fit of 
Eq. (HJ with A = 1.57 and ^jv = 0.35. Hence for {u) = 1.41 
we finds (3 = 0.49 in perfect agreement with the predicted 
value (3 = i/23 ,5,6,7,s,9,io,ii _ j^^^g^. rpj^^ critical disorder Wc 

as function of (j) for larger values of the magnetic fiux. The 
squares correspond to the values of Wc{4>) depicted in Tab.|I] 
and nil The red line is the power law fit used in in the main 
panel. = and (f> = 0.25 are symmetric axes^^ and therefore 
the whole range of (j> is covered. One can see the fast shift of 
W^c for < 3 X 10"^, which is the region for which the scaling 
holds. 



prediction of Larkin and Khmel'nitskii were verified it is 
nevertheless quite surprising how fast the critical disorder 
rises for small magnetic fiuxes. Already at = 0.03 most 
of the shift in the critical disorder has taken place. For 
higher magnetic fluxes (cj) > 0.03) the deviations from 
the power-law scaling are seen. We believe that such 
a sensitivity to small magnetic fluxes may be used as a 
basis for building a magnetic sensor. 



IV. THE CRITICAL LEVEL SPACING 
DISTRIBUTION 

In contrast to the critical disorder Wc the critical level- 
spacing distribution P{s) is not universal, and it depends 
on the boundary conditionsi^^ . In the absence of a mag- 
netic field^- the peak of the distribution is shifted to 
smaller values of s when the boundary conditions are 
successively changed from periodic to hard-wall in each 
of the three spatial directions. The tails of P{s) for large 
s can be fitted by an exponential P{s) ~ exp(— ks) with 
K — 1.87 for periodic boundary conditions— and k = 1.49 
for hard-wall (Dirichlet) boundary conditionsiSi. The 
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TABLE II: Results and errors of simultaneous fit for the interesting parameters. 



(0.0008) 


VKc(O.OOI) 


H^c(0.00I5) 


Wc(0.002) 


H^e (0.004) 


H^c(0.007) 


WciO.Ol) 




loc 


V 


best fit 16.97 


f7.03 


17.09 


17.15 


17.32 


17.45 


17.56 


17.97 


0.838 


1.41 


error O.IO 


O.II 


O.II 


0.11 


0.12 


0.11 


0.11 


0.12 


0.006 


0.04 




s 



FIG. 4: (Color online) The critical level spacing distribution 
functions for periodic and hard-wall boundary conditions and 
different fluxes </>. For periodic boundary conditions they de- 
pend on 0, where the dashed red line corresponds to (ji = 
and the dashed green line to = 0.1. For hard- wall bound- 
ary conditions, there is no such dependence, as can be seen 
from the behavior of the continuous lines from red to green 
corresponding to <;/> = 0, 8 x 10"", 10"^, 1.5 x 10"^, 2 x 10"^, 
4x 10-^ 7x 10-^ 10-^ 3x 10-^ O.l, and 0.25. The limiting 
distributions axe also shown, black dotted line - Poisson dis- 
tribution, black dashed line - GOE, black dashed dotted line 
- GUE. Inset: a semi-log plot depicting the behavior for large 
s. The solid blue lines are fits to the mean tails with k = 1.42 
in the case of hard-wall boundary conditions and k = 1.88 for 
periodic boundary conditions, both independent of 0. 

critical level spacing distribution does not only change 
with boundary conditions but also when a magnetic flux 
is applied. The shape shifts towards a more GUE-like 
shape for periodic boundary conditions^-. However, we 
are not aware of any numerical calculation for hard-wall 
boundary conditions. 

Figure [5] shows the critical level-spacing distribution 
functions P{s) for Anderson models with different mag- 
netic fluxes (j) and periodic as well as hard-wall bound- 
ary conditions. The results fully conflrm the previous 
studies discussed in the previous paragraph. However, 
the critical level spacing distribution P{s) is indepen- 
dent of (p. This new behavior of P{S) for hard-wall 
boundary conditions is surprising, since the universality 
class changes from 4> — (limit: GOE) to > (limit: 



GUE) as can be seen by the different v values (see Ta- 
ble I) . In addition, the critical level distribution function 
changes when periodic boundary conditions are applied. 
The tails of the critical distribution confirm to the ansatz 
P{s) exp(— Ks) with k = 1.88 for periodic boundary 
conditions and n = 1.42 for hard- wall boundary condi- 
tions independent of (j). The fits are depicted in the inset 
of Fig. H 

The critical second moment /o,c depends also on the 
boundary conditions. In the literature one finds values 
for periodic boundary conditions, /o.c = 0.7 1^^'^"^ and for 
hard-wall boundary conditions Iq^c = 0.81— for </> = 0. 
We confirm both numbers numerically for = and 
(/) 7^ 0. In the case of periodic boundary conditions the 
shape of P{s) changes when a magnetic flux is introduced 
but the critical second moment stays the same. 



V. SUMMARY AND CONCLUDING REMARKS 

Our main goal has been to confirm that the shift 
in the critical disorder Wc{(p) for small values of the 
magnetic flux (p follows the scaling behavior Wc{(j)) — 
Wc{0) = Acj)^^" with /3 = 1/2 as predicted by Larkin and 
Khmel'nitskii^. After a careful numerical study we are 
able to conflrm this prediction, however, the prefactor 
A is much larger than expected from previous numeri- 
cal studiesiS,. Thus, the scaling holds only for a narrow 
range of small fluxes up to (/> « 0.03 resulting in the 
fact that very weak fluxes have a strong effect on the re- 
sistivity of a sample at low temperature. This type of 
behavior naturally leads to the observation that the ef- 
fect might be used as the basis of an extremely sensitive 
low-temperature sensing device. 

In addition we have further considered the effects of 
boundary conditions and magnetic fields on the details 
of the metal-insulator transition of a 3D Anderson model. 
The results could be summarized as follows: the critical 
disorder Wc is independent of the boundary conditions 
although it depends on (j). The opposite is true for the 
critical second moment /q.c of the level spacing distribu- 
tion P{s), which is independent of (j) but changes with 
boundary conditions. The flux dependence of the critical 
level spacing P{s) is sensitive to the boundary conditions. 

Acknowledgment: This work has been supported by 
the Minerva Foundation, the Israel Science Foundation 
(Grant 569/07), and the Deutsche Forschungsgemein- 
schaft (DFG, within SFB 418). 



Handbook of Magnetism and Advanced Magnetic Materials, 
edited by H. Kronmiiller and S. Parkin, (John Wiley & 
Sons, 2007). 

P. W. Anderson, Phys. Rev. 109, 1492 (1958). 

D. E. Khmornitskii and A. I. Larkin, Solid State Commun. 

39, 1069 (1981). 

K. Slcvin and T. Ohtsuki, Phys. Rev. Lett. 78, 4083 (1997); 
T. Kawarabayashi, B. Kramer, and T. Ohtsuki, Phys. 
Rev. B 57, 11842 (1998); T. Ohtsuki, K. Slevin, and T. 
Kawarabayashi, Ann. Phys. (Leipzig) 8, 655 (1999). 
JY. Yuan and T. R. Kirkpatrick, J. Stat. Phys. 64, 309 
(1991). 

M. Biafore, C. Castellani, and G. Kothar, Nucl. Phys. B 
340, 617 (1990). 

C. Castellani and G. Kotliar, Physica A 167, 294 (1990). 
F. J. Wegner, Nucl. Phys. B 270, 1 (1996). 

R. Oppermann, J. Phys. Lett. (Paris) 45, 1161 (1984). 

D. Belitz, Solid State Commun. 52, 989 (1984). 
B. Shapiro, Philos. Mag. B 50, 241 (1984). 

T. Drose, M. Batsch, L K. Zharckcshcv, and B. Kramer, 
Phys. Rev. B 57, 38 (1998). 

V. F. Mitin, V. K. Dugaev, and G. G. Ihas, App. Phys. 
Lett. 91, 202107 (2007). 



M. Watanabc, K. M. Itoh, Y. Ootuka, and E. E. Haller, 
Phys. Rev. B 60, 15818 (1999). 

E. Tousson, V. Volterra, E. P. Rubenstein, R. Rosenbaum, 
and Z. Ovadyahu, Philos. Mag. B 56, 875 (1987). 

^® M. Batsch, L. Schweitzer, L Kh. Zharekeshev, and B. 
Kramer, Phys. Rev. Lett. 77, 1552 (1996). 

D. Braun, G. Montambaux, and M. Pascaud, Phys. Rev. 
Lett. 81, 1062 (1998). 

L. Schweitzer and H. Potempa, Physica A 266, 486 (1999). 
1" D. R. Hofstadter, Phys. Rev. B 14, 2239 (1976). 

B. L Shklovskii, B. Shapiro, B. R. Sears, P. Lambrianides, 
and H. B. Shore, Phys. Rev. B. 47, 11487 (1993). 

E. Hofstettcr and M. Sclircibcr, Phys. Rev. B 48, 16979 
(1993); ibidPhys. Rev. B 49, 14726 (1994). 

K. Slevin and T. Ohtsuki, Phys. Rev. Lett. 82, 382 (1999); 
K. Slevin, T. Ohtsuki, and T. Kawarabayashi, ibid84, 3915 
(2000). 

23 A. Stathopoulos, SIAM J. Sci. Comput. 29, 481 (2007); A. 

Stathopoulus and J. R. McCombs, ibid 29, 2162 (2007). 
2* L Varga, E. Hofstetter, M. Schreiber, and J. Pipek, Phys. 

Rev. B 52, 7783 (1995). 



